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Appareil a liquide pour l'intégration graphique de 
certains types d'équations différentielles. 

Par M. Michel Petrovitch h Belgrade {Serbie), 



Dans mon Article, " Sur l'intégration hydraulique des équations différen- 
tielles " (American Journ. of Mathem., vol. XX, No. 4), j'ai indiqué une méthode 
pouvant servir à intégrer graphiquement des équations différentielles du premier 
ordre. On peut faire varier la disposition de l'appareil et par cela même les 
types d'équations que la méthode permet d'intégrer. 

J'ai fait construire un appareil de ce genre, particulièrement facile à réaliser, 
pouvant servir à enregistrer directement les intégrales de certains types d'équa- 
tions. Ces équations sont, il est vrai, intégrables analytiquement, mais il est 
commode pour les applications, d'avoir une méthode rapide et sûre pour la con- 
struction mécanique de leurs courbes intégrales. 

Description et fonctionnement de Vappareil. 

Supposons que sur une plaque d'épaisseur uniforme (p. ex. sjir une planche 
en bois) on marque un système d'axes rectangulaires oz, ot (fig. 1) ; qu'on y trace 
ensuite une courbe donnée 

t-F{z) 
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et qu'on découpe la surface limitée par deux droites Zi et z^, la courbe ainsi tracée 
et l'axe des a (si la courbe s'étend de deux cotés de oa, on découperait la plaque 
suivant cette courbe et les deux droites z^ et z^. Cette partie découpée sera 
appelée dans la suite : plaque D. 

Trapons ensuite sur une autre plaque, d'épaisseur également uniforme, un 
système d'axes ov, ow et une courbe donnée 

et comme précédemment découpons la partie de la plaque comprise entre la 
courbe et les deux droites Vi, v^. Cette partie découpée sera appelée: 
plaque E. 

Supposons maintenant que sur un cylindre G on ait tracé deux axes 
0^ et ox, de sorte que o^ soit parallèle à la base du cylindre et que ox soit une 
des ses génératrices. Traçons y ensuite une courbe donnée 

x = e{^). 

Ceci étant, supposons un appareil construit d'après le schéma suivant. 

Une cuve ^(fig. 2), ayant la forme d'un parallelipipède, communique par 
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un tube horizontal M avec un autre tube vertical N. Le cylindre G est placé 
sur un support (p. ex. sur une caisse en bois) et peut tourner lentement autour 
d'un axe vertical soit par l'action d'un mécanisme d'horlogerie, soit par une 
manivelle qu'on tourne à la main. 

On fixe la plaque E de sorte qu'elle immerge dans la cuve ^et qu'elle ne 
puisse pas se déplacer pendant l'expérience. 

La plaque D, également immergée dans la cuve, est fixée au bout d'une tige 
P pouvant glisser verticalement. Sur la tige P se trouve fixé un repère B ; en 
le tenant par la main, on pourra, en le déplaçant verticalement ensemble avec la 
tige (et par suite ensemble avec la plaque D) faire de sorte que lorsque le 
cylindre G tourne, le repère se trouve constamment sur la courbe a; = (^) 
tracée d'avance sur le cylindre. 

Versons ensuite d'un liquide quelconque dans la cuve K, jusqu'à une hau- 
teur arbitraire. Lorsque le cylindre G commence à tourner et qu'en déplaçant 
à chaque instant le repère B on fait de sorte qu'il se trouve constamment sur la 
courbe a; = (^) , le niveau du liquide s'élèvera ou baissera avec la longueur ^ 
parcourue par un point du cylindre suivant une loi dépendant des formes 

1° de la courbe x=Q{^) tracée sur le cylindre ; 

2° de la courbe t = F{z), suivant laquelle a été découpée la plaque D ; 

3° de la courbe w=i^{v), suivante laquelle a été découpée la plaque E. 

Supposons qu'une tige verticale S soit placée à coté du tube N, pouvant 
glisser verticalement et sur laquelle se trouve fixe un repère T en bas et un 
crayon Ç à la hauteur du cylindre. Si un observateur déplace à chaque instant 
le repère l'de manière que sa pointe se trouve constamment à la hauteur du 
niveau du liquide dans le tube N et que, à l'aide d'une vis, on fixe le crayon Q 
sur la tige S à& sorte que la distance QT soït invariable pendant ce déplacement, 
le crayon décrira sur le cylindre tournant une certaine courbe 

dépendant des trois données 1°, 2°, 3° et dont nous allons chercher l'équation 
différentielle. 
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A cet effet désignons (fig. 3) par 
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a et 5 les épaisseurs uniformes des plaques D et E; 

Z et ^ la largeur et l'épaisseur de la cuve K; 

r le rayon du tube N; 

a la distance du repère B au point O, origine de coordonnées de la courbe 
suivant laquelle on a découpé la plaque D ; 

h la distance de l'origine de coordonnées pour la plaque E au plan du fond 
de la cuve K; 

/S la distance de l'axe /i/g des ^ sur le cylindre G au plan du fond de la 
cuve K-f 

p la distance du crayon Q au niveau du liquide ; 

S la distance de l'axe e^ e^ des ^ relatif à la courbe i/=zj^(^^) décrite par le 
crayon Q , au plan du fond de la cuve. On pourrait faire coïncider cet axe avec 
l'axe /1/2 relatif à la courbe décrite par le repère B. Mais comme ces deux 
courbes varient dans les sens inverses, il y a avantage dans la pratique, pour 
diminuer la hauteur du cylindre, de choisir convenablement la hauteur de cet axe ; 

X l'ordonnée de la courbe x=6{^) tracée sur la cylindre ; 

p l'ordonnée de la courbe y = j^(^) décrite par le crayon Q ; 

z la distance du point au niveau du liquide ; 
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u la hauteur de ce niveau comptée à partir du fond de la cuve K; 

V la distance de l'origine des coordonnées de la plaque E au niveau du 
liquide. 

Les quantités a, b, l, g,r, a, h, (3, p, è sont des constantes pour une expé- 
rience donnée ; les x, y, z, u, v sont des variables. 

Ceci étant, soit ES' le niveau initial du liquide. En faisant immerger la 
plaque D de sorte que x se change en a; — dx, le volume dV d\i liquide qui s'est 
élevé au-dessus du niveau primitif u sera 

dV = gl du, + r^ndu — b^ (v) du — ■ aF{z) d,u (1) 

et comme l'on a 

u = y + S-p, (2) 

v = h — M = A + p — ^ — 2/. (3) 

en posant pour abréger 

gl^^n-b^{h + p-8-y) = ^{y), (4) 

aF{z) = ^{z) (5) 

on aura 

dV=l^{y)-^{z)-]dy. (6) 

Mais ce volume est égal au volume du liquide déplacé par l'immersion de 
la plaque D. Comme la largeur de la partie immergée audessous du niveau 
HE' est F{z), son épaisseur a et la hauteur — dx, son volume sera 

— aF{z) dx = —^ (z) dx 
d'où l'équation 

[<|> (y) _ tp (z)] dy = — ^! (z) dx . (7) 

Mais on a à chaque instant 

u + z + a = x-^^ (8) 

d'où z = x — y + :?i,] /gx 

x^z + y — /l,i 

avec a=/3 + S — a — p; (10) 

par suite on peut écrire l'équation (7) sous l'une on l'autre des formes suivantes : 

^{y)dy + ^{z)dz=0, (11) 

l^{y}-'f{^-y + m^ + 'f{^-y + ^) = ^- (12) 
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Telle est l'équatioa différentielle du problème. Si 

y,{z,y,G) = (13) 

est l'intégrale générale de (11), celle de (12) sera 

>7(a; — 2/ + X,2/, (7)=0, (14) 

et inversement : si 

f{x,y, G) = 

est l'intégrale de (12), celle de (11) sera 

f{z + y-%,y, G) = 0. 

Par suite, si la courbe, tracée sur le cylindre et le long de laquelle on fait 
déplacer le repère B, est 

et si l'intégrale générale de (11) est (13) , la courbe 3/ = ;^ (0) tracée par le crayon 
Q sera 

Cl étant une valeur particulière de la constante d'intégration, déterminée p. ex. 
par la condition imposée à la courbe y = ;^ (^) de passer par un point donné 
^ = ^j, 2/ = 2/0 ; ceci donne les trois équations 

Xf>=^ih), Zo= «•0 — 2/0 + ^. >7(zo. 2/0. <^) = 0» 

déterminant la valeur de G. 

Remarquons aussi qu'étant données : 

1° la courbe ce = (^), tracée d'avance ; 

2° la courbe 2/ = ;c (0 tracée sur le cylindre par le crayon Q ; 

3° les valeurs initiales ^o. ^o, 2/0 (liées d'ailleurs par 

2o = 9 (^0) — yQ + ^ 

et dont deux, par conséquent, sont arbitraires), on déterminera la valeur de z 
correspondant à une valeur donnée de y (et inversement) de la manière suivante. 
Si y z=.y^ est la valeur donnée de y, la valeur correspondante de z est 
l'ordonnée d'un des points d'intersection de la courbe 

y = H^)-yi + ^ (15) 

avec les droites 
parallèles à l'axe des y. 
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Inversement, si % est une valeur donnée de s, la valeur correspondante de y 
est l'ordonnée d'un des points d'intersection de la courbe 

avec la courbe 

y = H^)—^i-\-^- (16) 

Les courbes (16) et (16) ne sont autres que la courbe a; = 0(^) déplacée 
dans la direction de l'axe des x. 

Application h quelques problèmes. 

I. D'abord, l'appareil peut servir comme intégrateur, permettant de calculer 
la valeur d'une intégrale définie 

J=f/{z)dz. 

A cet effet supposons ; 

1° qu'on a donné à la plaque E la forme d'un parallélipipède, de sorte 
qu'on ait 

<ï> (^) = const = B ; 

2° que la courbe x=6{^) soit une droite quelconque 

X ^ m^ + n ; 

S" qu'on a donné à la plaque D une forme telle, que la fonction précédente 
F (z) soit égale à /(s), c'est-à-dire qu'on ait 

^{z) = a/{z). 
L'intégrale de (11) sera alors 

z 
y=— ^J/(2)<^3+2/0 

et, d'autre part, la valeur yi, correspondant à la valeur z = %, sera l'ordonnée du 
point d'intersection de la droite 

y = m^ + %i 

(oùîli = n + ^ — Zi) avec la courbe ?/ = ;^(|), tracée par le crayon Q. Cette 
valeur y^ étant ainsi, déterminée, on aura 

«^= — —(2/1 — 2/0), 
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II. L'appareil peut servir aussi comme intêgraphe. A cet effet supposons 
1° qu'on a donné à la plaque D la forme d'un parallélipipède de sorte 

qu'on ait 

F{z) = const = H; 

2° que la courbe x = 0{^) soit une droite quelconque 

x=-m^ -\- n. 

L'équation différentielle (12) devient alors, après avoir posé 

dx = md^ , 
[$ {y) — aE] dy + maHd^ = 

d'oii 



/[CD (y) - aH2 dy=- maH{^ - ^„) . 



Par conséquent, d'après l'équation (4), pour que l'abscisse ^ de la courbe 
yz=.j^{^), tracée sur le cylindre par le crayon Q, représente à chaque instant la 
valeur d'une intégrale donnée 



J=^ff{y)dy, 



il faut et il suffit qu'on ait 



I/o 



aH—^(y) _ ... 



c'est-à-dire que la courbe, suivant laquelle a été découpée la plaque E, soit telle 
qu'on ait 

avec 

maH ql + t^ti — aH 7x1 s 

tti = — T— , «g = « ■ Y , « = ft + p . 

III. L'appareil s'applique à l'intégration graphique des équations 

f{y)dAj + -^{z)d^ = 0, (17) 

où /et '4' sont des fonctions positives données. 

Pour cela on donnera à la plaque D une forme telle qu'on ait 

i?'(,) = ±.^(.) 
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et à la plaque E une forme telle qu'on ait 

4>(2/) = A — /^s/C^ — 2/). 
avec 

enfin, on prendra pour le courbe œ = (^) une droite quelconque. D'après ce 
qui a été dit précédemment, en ayant la courbe ?/ = ;^ (^) tracée par le crayon Q 
et les valeurs initiales ^o, Zq, yo, on en déduit facilement la valeur y corres- 
pondant à une valeur donnée de a (ou inversement), y et z étant liées par l'équa- 
tion différentielle précédente. 

IV. L'appareil se prête aussi à l'intégration graphique des équations 

qui ne peuvent pas s'écrire directement sous la forme (17), mais dans lesquelles 
les variables peuvent être séparées par un changement de la variable indépen- 
dante de la forme 

^ = M(y + z), 

où z est la nouvelle variable indépendante et Jf une fonction donnée. 

En effet, en effectuant ce changement, on aura entre y et z une équation 

différentielle de la forme 

/{y)dy-{-^iz)dz=o 

(où nous supposerons les fonctions /et 4' positives). D'autre part, en désignant 

par N{E) la fonction inverse de 

^ = Mit), 

on aura N{^) = y -}- z = x + :>i = d{^) + ?^, 

d'où e{^) = N{^)-X, 

Par conséquent, pour que le crayon Q trace des courbes intégrales de (18), 
il faut et il suffit : 

1° qu'on donne à la plaque D une forme telle qu'on ait 

F{z) = ±^{z) 
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et à la plaque E une forme telle qu'on ait 

avec 

2° que la courbe se = ^ (^) soit la suivante : 

où /l a la signification précédente. 

V. L'appareil se prête à la résolution du problème suivant : étant donnée 

une fonction arbitraire 

t=:£l{t) 

et une courbe quelconque déjà tracée x= 6(^), tracer la courbe 

En remarquant que si cette courbe est celle décrite par le crayon Q, on aura 

et par suite y, définie comme fonction de z par cette équation, doit satisfaire à 
l'équation précédente 

<ï> {y) dy+^f (z) dz = 

(cette condition étant d'ailleurs suffisante). Cette condition ne dépend pas de la 
forme de la courbe x= d{^) et une fois qu'on l'a réalisée pour une quelconque 
des ces courbes, elle sera satisfaite pour une courbe quelconque. 

On la réalisera en donnant aux plaques D et U des formes convenables. 
Remarquons que la dernière condition revient à celle-ci : si 

t = A{r) 
est la fonction inverse de 

z définie par 

z=A{y)-y + 7^, 

doit satisfaire à l'équation (11), ce qu'on peut réaliser de bien de façons. P. ex. 
on donnera à la plaque D la forme d'un parallélipipède et à la plaque E une 
forme telle qu'on ait 

1-A'{y) = ^{y), 
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ou bien, en donnant à ^ la forme d'un parallélipipède de largeur 5 et à Z> une 
forme telle qu'on ait 

où l'équation y ^= j(^{z) représente la solution de 

y=a{z-\-y-'X) 
par rapport à y etc. 



Au lieu d'avoir une plaque mobile et une fixe, on peut en avoir plusieurs. 
Désignons les plaques mobiles par Z>j et les plaques fixes par Ei . Supposons 
que toutes les plaques Di se meuvent indépendemment les unes des autres et qu'à 
la plaque £>i correspond la courbe 

tracée sur le cylindre. En désignant comme précédemment par 

ti = Ft{z), 

la courbe suivant laquelle on ait découpé la plaque Di et par «j l'épaisseur de 
celle-ci; par 

Wi = ^i{Vi), 

la courbe suivant laquelle on a découpé Ut et par bi l'épaisseur de celle-ci, on 
trouve facilement que la variation du niveau du liquide, contenu dans la cuve K, 
lorsque le cylindre se met en mouvement, sera réglée par l'équation différen- 
tielle 

^{y)dy+'^%{zi)dZi=0, (19) 

où l'on a posé pour abréger 

* (y) = gl + 'r^7t—^ bi ^i {h — y) , 

%{z,) = aiFi{Zi), 

les hi étant certaines constantes. En désignant par 

^{y,zi,^ 2ft) = o, (20) 

la relation entre y et les Zi, exprimée par l'équation 

/<ï> {y) dy + ^/'P* (^i) dzi + const., (21) 
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la courbe décrite par le crayon Q est celle, dont l'équation 

s'obtient en éliminant z^, Sg, . . . . , % entre (20) et les h équations 

22 = ^2(6 — 2/ + ^2. l (22) 

2. = ^.(^) — y + ^*, ^ 

où "ki sont certaines constantes analogues à la constante % précédente. 

On peut ainsi, en partant d'une équation quelconque de la forme (19) et en 
y remplaçant les %i par leurs valeurs (22) (oii les fonctions 6» et les constants "k^ 
peuvent être quelconques), former des équations du premier ordre, en ^ et «/, 
dont les courbes intégrales peuvent être enregistrées par le crayon Q de 
l'appareil. 

Remarque. Je saisis l'occasion de signaler quelques erreurs de calcul, qui se 
sont glissé dans mon Article cité ; on rétablira aisément les formules exactes en 
remplaçant dans cet Article l'équation (l) par la suivante: 

l^{y)-F{z)\dy^-F{z)dx 

(comme il est facile de s'en rendre compte) et par suite en modifiant les formules 
qui en résultent. Cela, d'ailleurs, n'atteint pas la signification de la méthode et 
ne change que des types intégrables par les procédés décrits. 



